Chapter 2 

Second Quantisation 



In this section we introduce the method of second quantisation, the basic framework for 
the formulation of many-body quantum systems. The first part of the section focuses on 
methodology and notation, while the remainder is devoted to physically-motivated appli- 
cations. Examples of the operator formalism are taken from various fields of quantum 
condensed matter. 



2.1 Notations and Definitions 

Second quantisation provides a basic and efficient language in which to formulate many- 
particle systems. As such, extensive introductions to the concept can be found throughout 
the literature (see, e.g., Feynman's text on Statistical Mechanics [10]). The first part 
of this section will be concerned with the introduction of the basic elements of second 
quantisation, while the remainder of this section will be concerned with developing fluency 
in the method by addressing a number of physical applications. 

Let us begin by defining the (normalised) wavefunctions and corresponding eigen- 
values e\ of the single-particle Hamiltonian H, viz. 

H\1>\) =ca|^a>- 

With this definition, populating states 1 and 2, the symmetrised (normalised) two-particle 
wavefunction for fermions and bosons is respectively given by 

i/jf(x 1 ,x 2 ) = -j= (^i(x 1 )i) 2 {x 2 ) T ^2(^1)^1(^2)) • 



In the Dirac bracket representation, we can write 

\1,2)f = -±= (|V>i> ® T 1^2) ® ■ 

More generally, a symmetrised iV-particle wavefunction of fermions (( = —1) or bosons 
(C = +1) is expressed in the form 
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where n\ is the total number of particles in state A (for fermions, Pauli exclusion en- 
forces the constraint n\ = 0, 1, i.e. ri\\ = 1) - see Fig. 2.1. The summation runs over 
all iV! permutations of the set of quantum numbers {Ai,...Ajv} 5 and V denotes the 
parity, defined as the number of trans- 
positions of two elements which brings 
the permutation (Vi,V2, ■ • • Vn) back 
to the ordered sequence (1, 2, • • • N). 
Note that the summation over per- 
mutations is necessitated by quantum 
mechanical indistinguishability: for 
bosons/fermions the wavefunction has 

to be symmetric/anti-symmetric under particle exchange. It is straightforward to confirm 



Enrico Fermi 1901-1954: 1938 
Nobel Laureate in Physics for his 
demonstrations of the existence of 
new radioactive elements produced 
by neutron irradiation, and for his 
related discovery of nuclear reactions 
brought about by slow neutrons. 




that the prefactor 



normalises the many-body wavefunction. In the fermionic 



case, the many-body wavefunction is known as a Slater determinant. 
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Figure 2.1: Schematic showing typical occupation numbers for a generic fermionic and bosonic 
system. 

The expression above makes it clear that this 'first quantised' representation of the 
many-body wavefunction is clumsy. We will see that the second quantisation provides the 
means to heavily condense the representation. Let us define the vacuum state and 
introduce a set of field operators a\ together with their adjoints at, as follows: 1 



a A |0> = 0, 



a Ail^) = |Al,A 2 ,...Ajv) 



(2.1) 



Physically, the operator a\ creates a particle in state A while the operate a\ annhilates it. 
These definitions are far from innocent and deserve some qualification. Firstly, in order 
not to be at conflict with the symmetry of the wavefunction, the operators a\ have to 
fulfill the commutation relations, 



— 5\,u, 



= 0, 



-c 



= 



-c 



(2.2) 



x As before, it will be convenient to represent these operators without a circumflex. 
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where [A, B]_^ = AB—(BA is the commutator ( = 1 (anticommutator ( = — 1) for bosons 
(fermions). 2 The most straightforward way to understand this condition is to check that 
the definition |A, /u) = a\a^\fl) and property |A,/i) = C,\/i, A) in fact necessitate Eqs. (2.2). 
Yet even if (2.2) is understood, the definitions above remain non-trivial. Actually, quite a 
strong statement has been made: for any N, the A-body wavefunction can be generated 
by an application of a set of N -independent operators to a unique vacuum state. In order to 
check that Eqs. (2.1) and (2.2) actually represent a valid definition, including, for instance, 
the right symmetrisation and normalisation properties of A^-body wave functions, various 
consistency checks have to be made. 

Based on Eqs. (2.1) and (2.2), a formal definition of the general many-body or Fock 
space can now be given as follows. First define Tn to be the linear span of all N- 
particle states |Ai, ■ • ■ A^r) = —j= a\ N ■ ■ ■ a^JQ). The Fock space T is then defined as 

the direct sum ©^=0-^ ( see Fig- 2. 2). 3 A general state \4>) of the Fock space is, therefore, 




Figure 2.2: Visualisation of the generation of the Fock-subspaces Fn by repeated action of 
creation operators onto the vacuum space J-q. 



a linear combination of states with any number of particles. To turn these rather abstract 
definitions into a valuable tool for practical computation we need to put them into relation 
with standard operations performed in quantum mechanics. In particular we have to 
specify how changes from one single-particle basis {A} to another {A} effect the operator 
algebra {a\}, and in what way standard operators of (many-body) quantum mechanics 
can be represented in terms of the a\ s: 



> Change of basis: Using the resolution of identity, id = YlT=o l^)(^l> the relations |A) = 
| A) (A | A) , | A) = a\\Q,), and |A) = a~|f2) immediately give rise to the transformation 

law 



4 = EA<A|A)a 



Ea(A|AK 



(2.3) 



In many applications we are not dealing with a set of discrete quantum numbers (spin, 
quantised momenta, etc.), but rather with a continuum (a continuous position coordinate, 
say). In these cases, the quantum numbers are commonly denoted in a bracket notation 
a\ ~> a(x) = ^^ A (x| A)aA, and the summations appearing in the transformation formula 
above become integrals. 



2 As a convention, when unspecified by the notation [•, •] will be used to denote the commutator and 
{•, •} the anticommutator. 

Here, the symbol of the direct sum © is used to show that each "submodule" Tn is linearly indepen- 
dent. 
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Example: The transformation from the coordinate to the momentum representation 
in a finite one-dimensional system of length L would read 



dk = dx (k\x)a(x), a(x) = 'S^(x\k)a,k, 
Jo u 



where (k\x) = (x\k)* = jjj^e , cf. Fourier series expansion. 

> Representation of operators (one-body): Single particle or one-body operators 
Ox acting in a iV-particle Hilbert space, Tn, generally take the form 0\ = ^2„ =1 6 n , 
where d n is an ordinary single-particle operator acting on the n-th particle. A typical 
example is the kinetic energy op- 



erator T = ^ n ^, where p n is 
the momentum operator acting on 
the n-th particle. Other examples 
include the one-particle potential 
operator V = ^2 n V(x n ), where 
V(x) is a scalar potential, the 



David Hilbert 1862-1943: His work in 
geometry had the greatest influence in 
that area after Euclid. A systematic study 
of the axioms of Euclidean geometry led 
Hilbert to propose 21 such axioms and 
he analysed their significance. He con- 
tributed to many areas of mathematics. 




total spin-operator ^ n S„, etc. 

Since we have seen that, by applying field operators to the vacuum space, we can gener- 
ate the Fock space in general and any iV-particle Hilbert space in particular, it must be 
possible to represent any operator (D\ in an a-representation. 

Now, although the representation of n-body operators is after all quite straightforward, 
the construction can, at first sight, seem daunting. A convenient way of finding such a 
representation is to express the operator in terms of a basis in which it is diagonal, and 
only later transform to an arbitrary basis. For this purpose it is useful to define the 
occupation number operator 



n\ = a[a x 



(2.4) 



with the property that, for bosons or fermions (exercise), fix (atj"|Q) = n (a\) n \Q,), i.e. 
the state (a\) n \Q,) is an eigenstate of the number operator with eigenvalue n. When acting 
upon a state |Ai, A2, • • • Ajv), it is a straightforward exercise to confirm that the number 
operator simply counts the number of particles in state A, 



h\\\i, A 2 , • • • Xn) = a\a\ 



N 



' a Ail^) = 5^^AA i |Ai,A 2 ,---A iV ) 



i=l 



Let us now consider a one-body operator, 0\, which is diagonal in the orthonormal 
basis |A), 6 = Xa°a|A)(A|, o\ = (A|o|A). With this definition, one finds 



N 



(A'd • • • X' N \Oi\\i, ■ ■ ■ \ N ) — o\i j (A'j, • • • A'jyjAi, • • • X N ] 

= (X[, ■ ■ ■ X' N \y^oxnx\X 1: ■ ■ ■ Ajv). 



A=0 
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Since this equality holds for any set of states, we obtain the operator or second quantised 
representation 0\ = YlxLo ^ x = ^2T=o(^\^) a \ a ^- The result is straightforward; a 
one-body operator engages a single particle 



at a time — the others are just spectators. 
In the diagonal representation, one simply 
counts the number of particles in a state A 
and multiplies by the corresponding eigen- 
value of the one-body operator. Finally, by 
transforming from the diagonal representa- 
tion to a general basis, one obtains the result, 



Wolfgang Pauli 1900-1958: 
1945 Nobel Laureate in Physics 
for the discovery of the Exclu- 
sion Principle, also called the 
Pauli Principle. 




6i = 5>|A> OA <A|i/>ata,, 









(2.5) 



Formally, the one-body operator, Oi, scatters a particle from a state v into a state \i with 
probability amplitude (n\5\v). 

Examples: The total spin operator is given by 



aa ' a Xa' ' 



Acta' 



(2.6) 



where a =t, 4- is the spin quantum number, A denotes the set of additional quantum 
numbers (e.g. coordinate), and a denotes the vector of Pauli spin matrices 



<Jt. = 



1 

1 



<7„ 



-i 

1 



1 

-1 ' ' 



(2.7) 



i.e. S z = \ ^ A (n At - n X i), and 5+ = ^ A 4t a A|- 

Second quantised in the position representation, the one-body Hamiltonian is given 
as a sum of kinetic and potential energy as (exercise) 



H = T + V 



-I 



dx a\a 



P 

2m 



V(x) 



a(x) 



where p = —ihd x . (Note that the latter is easily proved by expressing the kinetic energy 
in the diagonal (i.e. momentum) representation — see problem set.) 

Finally, the total occupation number operator is defined as N = J dx a*(x)a(x). 

> Representation of operators (two-body): Two-body operators (D 2 are needed to 
describe pairwise interactions between particles. Although pair-interaction potentials are 
straightforwardly included into classical many-body theories, their embedding into con- 
ventional many-body quantum mechanics is made awkward by particle indistinguishabil- 
ity. As compared to the conventional description, the formulation of interaction processes 
within the language of second quantisation is considerably more straightforward. 

Initially, let us consider particles subject to the symmetric two-body potential V(x, x') = 
V(x',x). Acting on two-particle states, the operator is given by 



y(2) 



dx' \x,x')V(x,x')(x,x'\. 



(2i 
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Our aim is to find an operator V in second quantised form whose action on a many-body 
state gives 



TV ^ TV 

V\x 1 ,x 2 , ■ ■ ■ x N ) = V(x n ,x m )\x 1 ,x 2 , ■ ■ ■ x N ) = o ^ V(x n ,x m )\xi,x 2 , ■ --xn). 

Comparing this expressions with (2.8) one might immediately guess that 

V = ^JdxJdx' a\x)a i (x')V{x,x')a(x')a{x). 

That this is the correct answer can be confirmed by applying the operator to a many-body 
state. We first note that 

TV 

a(x')a(x)\x 1 ,x 2 , ■ ■ -x N ) = a(x')'y]C~ 1 ^(x - x n )\x 1 ,x 2 , • • • (no x n ) ■■■x N ) 

n=l 

TV TV 

= ^2C~ l S(x - x n ) ^2 VmnS(x' - x m )\x 1 ,x 2 , ■ ■ ■ (no X n ,X m ) ■ ■ -X N ) 

n=l m=l,(ro/ri) 



where 



( m 1 if ;// • ;/ 
[ Q if m > n 

Then, making use of the identity 

a< (x)a^ (x')a(x')a(x)\xi, x 2 , ■ ■ ■ x^) 

TV 

= ^2 C^VmJix - x n )5(x - x m )\x,x',x 1 ,x 2 , • • • (no x n ,x m ) ■■■xn) 

TV 

^ ^ C Vmn^^X X n )(5(x X m ) \x n , X m , X\ , X 2 , (no X ni X m ) • • • Xtv) 
TV 

= y^S(x- x n )S(x' - x m )\x 1 ,x 2 , ■ ■ -x N ), 

m^=n 

multiplying by V(x, x')/2, and integrating over x and x' , one confirms the validity of the 
expression. It is left as an exercise to confirm that the expression, | J dx J dx' V(x, x')n{x)n{x') 
although a plausible candidate, does not reproduce the two-body operator. 

More generally, turning to a non-diagonal basis, it is easy to confirm that a general 
two-body operator can be expressed in the form 



62= °w',w a \> a W a \> ( 2 - 9 ) 
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where 0^ xx > = (fJ., fjf\0 2 \X, X'). 

In principle one may proceed in the same manner and represent general n-body in- 
teractions in terms of second quantised operators. However, as n > 2 interactions rarely 
appear, we refer to the literature for discussion. 

This completes our formal introduction to the method of second quantisation. To 
make these concepts seem less abstract, the remainder of this section is concerned with 
the application of this method to a variety of problems. 

2.2 Applications of Second Quantisation 

Although the second quantisation is a representation and not a solution, its application 
often leads to a considerable simplification of the analysis of many-particle systems. To 
emphasize this fact, and to practice the manipulation of second quantised operators, we 
turn to several applications. The first example is taken from the physics of correlated 
electron systems, and will engage the manipulation of fermionic creation and annihilation 
operators. The second example involves the study of quantum magnetism within the 
framework of boson creation and annihilation operators. However, before getting to these 
applications, let us first go back and reinterpret our analysis of phonon modes in the 
quantum chain. 

2.2.1 Phonons 

Although, at the time, we did not specify in which Hilbert space the field operators 
act, the answer is that the representation space is again a Fock space; this time a Fock 
space of phonons or, more formally, of oscillator states. In contrast to what we'll find for 
the fermion case below, the Fock space in the phonon problem does not have an a priori 
interpretation as a unification of physical A-particle spaces. However, outgoing from 
a vacuum state, it can be constructively generated by applying the oscillator creation 
operators a\ to a unique vacuum state: 

> INFO. Define a ground or vacuum state by requiring that all operators annihilate 
it. Next define to be the space generated by We may then introduce a set of states 
\k) = aJjQ), k = 0, 2n/L, ... by applying oscillator creation operators to the vacuum. Physically, 
the state \k) has the significance of a single harmonic oscillator quantum excited in mode k. In 
other words, all oscillator states k! ^ k are in their ground state, whilst mode k is in the first 
exited state. The vector space generated by linear combinations of states \k) is called T\. This 
procedure can be iterated in an obvious manner. Simply define the space Tn to be generated by 
all states . . . a\. |0) = \k\, . . . , kw)- The spaces can be defined more concisely by saying 
that they are the eigenspaces of the occupation number operator with eigenvalue N. Finally, 
the Fock space is just the direct sum of all J-n, J- = ©^ = gJ r 7v- By construction, the application 
of any one a' k or afc to states £ J- does not leave T . A closer analysis actually shows that the 
corresponding Fock space T represents a proper representation space for the operators afc. A 
particle interpretation of the phonon states can now be naturally introduced by saying that the 
Fock space sector Tn represents a space of bosonic A^-particle states. Application of a\ (a k ) to 
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a state 6 J-n creates (annihilates) a particle (cf. Fig. 2.2). 



2.2.2 Interacting Electron Gas 

As a second example, we will cast the Hamiltonian of the interacting electron gas in 
second quantised form. To emphasize the utility of this approach, in the next section we 
will use it to explore the phase diagram of a strongly interacting electron gas. In doing 
so, we will uncover the limitations of the "nearly free electron theory" of metals. 

As we have seen, in second quantised notation, the non-interacting Hamiltonian of a 
one-dimensional system of electrons subject to a lattice potential is given by 

where the fermionic electron field operators obey the anticommutation relations [c a (x), <? a i{x' 
8 [x — x') 8 aa i . The field operators act on the 'big' many-particle Fock space, T = 
®n=o^n- Each iV-particle space Tn is spanned by states of the form cj. (arjv) • ■ • cl i (xi)\fl) 
where the 'no-particle' state or vacuum |f2) is annihilated by all operators c CT (x). 

Applying a two-body Coulomb interaction potential, \ Y^i^j \ x C -x-\ > wnere x i denotes 
the position of the i-th electron, the total many-body Hamiltonian takes the second quan- 
tised form 





^ dx 


f dx ^4(x)4,(x') 

era' 


e 2 


-cA x ) c *( x ) 



> EXERCISE. Setting V(x) = and switching to the Fourier basis, reexpress the Coulomb 



2m 



+ V(x) 



c a (x), 



interaction. Show that the latter is non-diagonal, and scatters electrons between different quasi- 
momentum states — see Fig. 2.3. 

Having introduced both the field operators 
themselves and their representation spaces, we are 
in a position to point out certain conceptual analo- 
gies between the model theories discussed above. u' n-^V A „ 
In each case we have described a physical system \ v \ ( l) / 
in terms of a theory involving a continuum of op- Xy\y\y\. 
erators, <j>{x) (phonons) and c a (x) (electrons). Of k'+q,a'^ Ttk-q,a 
course there are also important differences between ' \ 
these examples. Obviously, in the phonon theory, 

we are dealing with bosons whilst the electron gas Figure 2.3: Feynman diagram- 
is fermionic. However, by far the most important matic representation of the two-body 
difference is that the first example has been a free Coulomb interaction, 
field theory. That means that the Hamiltonian 

contained field operators at quadratic order but no higher. As a rule, free field theo- 
ries can be solved (in a sense that will become clear later on) straightforwardly. The 
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fermionic model, however, represents a typical example of an interacting field theory. 
There are terms of fourth order in the field operators which arose from the Coulomb 
interaction term. As compared to free theories, the analysis of interacting theories is 
infinitely harder, a fact that will surely become evident later on. 

To develop some fluency in the manipulation of second quantised field operators we 
will continue by exploring the 'atomic limit' of a strongly interacting electron gas. In 
doing so, we will derive a model Hamiltonian which has served as a paradigm for the 
study of correlated electron systems. 

2.2.3 Tight-binding theory and the Mott transition 

According to the conventional band picture of non-interacting electrons, a system with 
a half-filled band of valence electron states is metallic. However, the strong Coulomb 
interaction of electrons can induce a phase transition to a (magnetic) insulating electron 
'solid' phase (much as interactions can drive the condensation of a classical liquid into a 
solid). To explore the nature of this phenomenon, known as the Mott transition after 
Sir Neville Mott (formerly of the Cavendish Laboratory), it is convenient to reexpress the 
interacting Hamiltonian in a tight-binding approximation. 




Figure 2.4: Infinitely separated, each lattice site is associated with a set of states, s = 0, 1, ■ ■ ■, 
bound to the ion core. Bringing together just two atoms, the orbitals weakly overlap and hy- 
bridise into bonding and anti-bonding combinations. Bringing together a well-separated lattice 
of atoms, each atomic orbital broadens into a delocalised band of Bloch states indexed by a 
quasi-momentum k from the Brillouin zone and an orbital or band index s. 

To develop an effective Hamiltonian of the strongly interacting electron system we 
begin by considering a lattice of very widely spaced (almost isolated) atoms — the atomic 
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limit. In the simplest non-interacting picture, the overlap of the outermost electrons 
(albeit exponentially weak) leads to a hybridisation of the electronic orbitals and leads 
to the 'delocalisation' of a narrow band of extended states (see Fig. 2.4). These Bloch 
states ipks{%) of H^°\ which carry a 

Felix Bloch 1905-1983: 1952 No- 
bel Laureate in Physics for the de- 
velopment (with Edward M. Purcell) 
of new methods for nuclear magnetic 
precision measurements and discov- 
eries in connection therewith. 



quasi-momentum index k and a band 
or orbital index s — 0, 1, • • •, provide a 
convenient basis with which to expand 
the interaction. We can, in turn, de- 
fine a set of local Wannier orbitals 




B.Z. 

£ 

kE[— 7r/a,7r/a] 



Akna 




\lpks) 



1 N 



-ikna 



\lpns) 



where N denotes the total number of primitive lattice sites (with periodic boundary 
conditions). Here the sum on k runs over the N Appoints spanning the Brillouin zone, i.e. 
k = 2irm/aN with integers —N/2 < m < N/2. 

If the lattice is very widely spaced, the 
Wannier state ip ns will differ little from the 
s-th bound state of an isolated atom at 
x = na (see Fig. 2.5). Restricting attention 
to the lowest band s = 0, and restoring the 
spin degrees of freedom a, the field opera- 
tors associated with the Wannier functions 
are defined by 




fn-Da 



(n+na 



L=Na 



dx 



cl(x)\n) 

\x) (x\lp n ), 



Figure 2.5: Diagram illustrating the weak over- 
lap of Wannier states in the atomic limit. 



i.e. 



N 



dx ip n (x)cl(x), 4(x) = ^2 V£( 



x CL 



71=1 



(2.11) 



Physically, can be interpreted as an operator creating an electron with spin a at site 
n in the lowest band. Since the transformation (2.11) is unitary, it is straightforward 
to confirm that the operators c ncr , and $ n(J obey fermionic anticommutation relations 



C mcr'\ + 



5 n 



In the presence of a two-body Coulomb interaction, a substitution of the field operators 
in Eq. (2.10) by Wannier operators generates the generalised tight-binding Hamiltonian 

(exercise) 



H 



4 



mnrs <ja' 



C mcr C n<j' C rcr' C scr' 



,,t 
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where the "hopping" matrix elements are given by 

tmn = -^ m \H {0) \^n) = £ ^'^k = C, 

k 

and the interaction parameters are set by (exercise) 

1 f L f L e 2 

Umnrs = " / cfo / rfx' ll)* m {x)lpl(x' ) Vv (x') 1p s (x) . 

^ Jo Jo \x — X \ 

Physically t mn represents the probability amplitude for an electron to transfer (hop) from 
a site m to a site n. 



> EXERCISE. Show that, in the Fourier basis, c^ a = ^= J2 n e mfc0£ W> * ne non-interacting 
Hamiltonian takes the diagonal form = Z e kc\ a Ck a - 

Expressed in the Wannier basis, the representation above is exact (at least for states 
contained entirely within the lowest band). However, for a widely spaced lattice, most of 
the matrix elements of the general tight-binding model are small and can be neglected. 
Focusing on the most relevant: 

> The direct terms U mnnm = V mn involve integrals over square moduli of Wannier 
functions and couple density fluctuations at different sites, 

where n m = Y^<j c ma c m(r- Such terms have the capacity to induce charge density 
instabilities. Here we will focus on transitions to a magnetic phase where such 
contributions are inconsequential and can be safely neglected. 

> A second important contribution derives from the exchange coupling which in- 
duces magnetic spin correlations. Setting J^ n = U mnmn , and making use of the 
identity a a p ■ a 7 s = 25^5^ — Sa^S^s, one obtains (exercise) 

^ ] ^ ] UmnmnCino C n<j' C m(T' C n<j = ~ ^ ^ ] ^ mn ( ' ^rt '^ l mJ l n J • 

m^n a a' m^n ^ ' 

Such contributions tend to induce weak ferromagnetic coupling of neighbouring 
spins (i.e. J F > 0). Physically, the origin of the coupling is easily understood as 
deriving from a competition between kinetic and potential energies. By aligning 
with each other and forming a symmetric spin state, two electrons can reduce their 
potential energy arising from their mutual Coulomb repulsion. To enforce the anti- 
symmetry of the two-electron state, the orbital wavefunction would have to vanish 
at x = x' where the Coulomb potential is largest. This mechanism is familiar from 
atomic physics where it is manifest as Hund's rule. 



Quantum Condensed Matter Field Theory 



30 



CHAPTER 2. SECOND QUANTISATION 



> However, in the atomic limit where the atoms are well-separated and the overlap 
between neighbouring orbitals is weak, the matrix elements Uj and Jf- are expo- 
nentially small in the interatomic separation. By contrast the 'on-site' Coulomb or 
Hubbard interaction 

y ] y ] ^rnmmm c l n(7 C ln(T' C ma' C m<7 = U ^mf^rrei; 
ma m 

where U = 2U mmmm , increases as the atomic wavefunctions become more localised. 

Therefore, dropping the constant energy off-set e = t nn , in the atomic limit, a strongly 
interacting many-body system of electrons can be described effectively by the (single- 
band) Hubbard Hamiltonian 



{mn) a m 



(2.12) 



where we have introduced the notation (mn) to indicate a sum over neighbouring lattice 
sites, and t = t mn (assumed real and usually positive). In hindsight, a model of this struc- 
ture could have been guessed on phenomenological grounds from the outset. Electrons 
tunnel between atomic orbitals localised on individual lattice sites and experience a local 
Coulomb interaction with other electrons. 

Deceptive in its simplicity, the Hubbard model is acknowledged as a paradigm of strong 
electron correlation in condensed matter. Yet, after forty years of intense investigation, 
the properties of this seemingly simple model system — the character of the ground state 
and nature of the quasi-particle excitations — is still the subject of heated controversey 
(at least in dimensions higher than one — see below). Nevertheless, given the importance 
attached to this system, we will close this section with a brief discussion of the remarkable 
phenomenology that is believed to characterise the Hubbard system. 

As well as dimensionality, the phase behaviour of the Hubbard Hamiltonian is char- 
acterised by three dimensionless parameters; the ratio of the Coulomb interaction scale 
to the bandwidth U/t, the particle density or filling fraction n (i.e. the average number 
of electrons per site), and the (dimensionless) temperature, T/t. The symmetry of the 
Hamiltonian under particle-hole interchange allows one to limit consideration to densities 
in the range < n < 1 while densities 1 < n < 2 can be inferred by 'reflection'. 

Focussing first on the low temperature system, in the dilute limit n < 1, the typical 
electron wavelength is greatly in excess of the particle separation and the dynamics is free. 
Here the local interaction presents only a weak perturbation and one can expect the prop- 
erties of the Hubbard system to mirror those of the weakly interacting nearly free electron 
system. While the interaction remains weak one expects a metallic behaviour to persist. 
By contrast, let us consider the half-filled system where the average site occupancy is 
unity. Here, if the interaction is weak U/t <^ 1, one may again expect properties remi- 
niscent of the weakly interacting electron system. 4 If, on the other hand, the interaction 
is very strong U/t^> 1, site double occupancy is inhibited and electrons in the half-filled 



4 In fact, one has to exercise some caution since the commensurability of the Fermi wavelength with 
the lattice can initiate a transition to an insulating spin density wave state characterised by a small 
quasi-particle energy gap — the Slater Instability 
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Figure 2.6: Conductivity of Cr-doped 
V2O3 as a function of decreasing pres- 
sure and temperature. At tempera- 
tures below the Mott-Hubbard transi- 
tion point (P c = 3738bar,T c = 457.5K) 
the conductivity reveals hysteretic be- 
haviour characteristic of a first order 
transition. Reproduced from Limelette 
et al, Universality and critical behavior 
at the Mott transition, Science 302, 89 
(2003). 



system become 'jammed': migration 
of an electron to a neighbouring lattice 
site necessitates site double occupancy 
incurring an energy cost U. In this 
strongly correlated phase, the mutual 
Coulomb interaction between the elec- 
trons drives the system from a metal 
to an insulator. 



Sir Neville Mott 1905-1996: 1977 
Nobel Laureate in Physics (with 
Philip W. Anderson and John H. van 
Vleck) for their fundamental theoret- 
ical investigations of the electronic 
structure of magnetic and disordered 
systems. 




t> INFO. Despite the ubiquity of the experimental phenomenon (first predicted in a cele- 
brated work by Mott) the nature of the Mott— Hubbard transtion from the metallic to the 
insulating phase in the half-filled system has been the subject of considerable debate. In the 
original formulation, following a suggestion of Rudolf Peierls, Mott conceived of an insulator 
characterised by two 'Hubbard bands' with a bandwidth ~ t separated by a charge gap U. 5 
States of the upper band engage site double occupancy while those states that make up the 
lower band do not. The transition between the metallic and insulating phase was predicted to 
occur when the interaction was sufficiently strong that a charge gap develops between the bands. 
Later, starting from the weakly interacting Fermi-liquid, Brinkman and Rice 6 proposed that the 
transition was associated with the localisation of quasi-particles created by an interaction-driven 
renormalisation of the effective mass. Finally, a third school considers the transition to the Mott 
insulating phase as inexorably linked to the development of magnetic correlations in the weak 
coupling system — the Slater instability. 



To summarise, we have shown how the method of second quantisation provides a 
useful and efficient way of formulating and investigating interacting electron systems. In 
the next section we will employ methods of second quantisation involving bosonic degrees 
of freedom to explore the collective excitations of quantum magnets. 

5 N. F. Mott, Proc. Roy. Soc. A 62, 416 (1949) — for a review see, e.g. N. F. Mott, Metal-Insulator 
transition, Rev. Mod. Phys. 40, 677 (1968) or N. F. Mott, Metal-Insulator Transitions, 2nd ed. (Taylor 
and Francis, London, 1990). 

6 W. Brinkman and T. M. Rice, Application of Gutzwiller's variational method to the metal-insulator 
transition, Phys. Rev. B 2, 4302 (1970). 
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2.2.4 Quantum Spin Chains 

In the previous section, emphasis was placed on charging effects generated by Coulomb 
interaction. However, as we have seen, Coulomb interaction may also lead to the indirect 
generation of magnetic interactions of both ferromagnetic and antiferromagnetic charac- 
ter. To address the phenomena brought about by quantum magnetic correlations, it is 
instructive to begin by considering systems where the charge degrees of freedom are frozen 
and only spin excitations remain. Such systems are realized, for example, in Mott insu- 
lators where magnetic interactions between the local moments of localised electrons are 
mediated by virtual exchange processes between neighbouring electrons. Here, one can 
describe the magnetic correlations through models of localised quantum spins embedded 
on lattices. We begin our discussion with the ferromagnetic spin chain. 



Quantum Ferromagnet 

The quantum ferromagnetic chain is specified by the Heisenberg Hamiltonian 



H = 



-j y s m • s 



m+l 



(2.13) 



where J > 0, and S m represents the quantum mechanical spin operator at lattice site m. 
In section 2.1 (cf. Eq. (2.6)) the quantum 



Werner Heisenberg 1901- 
1976: 1932 Nobel Laureate in 
Physics "for the creation of 
quantum mechanics, the appli- 
cation of which has, inter alia, 
led to the discovery of the al- 



mechanical spin was represented through 
an electron basis. However, one can con- 
ceive of situations where the spin sitting at 
site m is carried by a different object (e.g. 
an atom with non-vanishing magnetic mo- 
ment) . At any rate, for the purposes of our lotropic forms of hydrogen" . 
present discussion, we need not specify the 

microscopic origin of the spin. All we need to know is (i) that the lattice operators 5* 
obey the SU(2) commutator algebra (for clarity, we have set h~ = 1 in this section) 




(2.14) 



characteristic of quantum mechanical spins, and (ii) that the total spin at each lattice site 
is S. 7 

Now, due to the positivity of the coupling constant J, the Hamiltonian favours con- 
figurations where the spins at neighbouring sites are aligned in the same direction (cf. 
Fig. 2.7). A ground state of the system is given by |Q) = (8> m |<S m ), where \S m ) represents 
a state with maximal spin-z component: S^\S m ) = S\S m ). We have written 'a' ground 
state instead of 'the' ground state because the system is highly degenerate: A simulta- 
neous change of the orientation of all spins does not change the ground state energy, i.e. 
the system posesses a global spin rotation symmetry. 



7 Remember that the finite-dimensional representations of the spin operator are of dimension 25+1 
where S may be integer or half integer. While a single electron has spin 5 = 1/2, the total magnetic 
moment of electrons bound to an atom may be much larger. 
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Figure 2.7: Schematic showing the spin con- 
figuration of an elementary spin-wave excita- 
tion from the spin polarized ground state. 



> EXERCISE. Compute the energy expectation value of the state Defining global 

spin operator as S = ^ m S m , consider the state \a) = exp(i(n/2)a ■ S)|0). Making use of 

the Baker-Hausdorff identity, e ixd Ae~ ixd = A + iX[6,A] + ^[6, [6, A}} H or otherwise, 

verify that the state \a) is degenerate with Explicitly compute the state | (1, 0, 0)) . Con- 
vince yourself that for general a, \a) can be interpreted as a state with rotated quantisation axis. 

As with our previous examples, we expect that a global continuous symmetry will 
involve the presence of energetically low-lying exciations. Indeed, it is obvious that in 
the limit of long wavelength A, a weak distortion of a ground state configuration (cf. 
Fig. 2.7) will cost vanishingly small energy. To quantitatively explore the physics of these 
spin— waves, we adopt a 'semi-classical' picture, where the spin S 3> 1 is assumed to 
be large. In this limit, the rotation of the spins around the ground state configuration 
becomes similar to the rotation of a classical magnetic moment. 

[> INFO. To better understand the mechanism behind the semi— classical approximation, 

consider the Heisenberg uncertainty relation, AS l AS : > < | *S'- 7 ]) | = e* jfe |(5 fc )|, where AS' is 
the root mean square of the quantum uncertainty of spin component i. Using the fact that 
\(S k )\ < S, we obtain for the relative uncertainty, AS^S, ^f^f- < Jr ^> 0, i.e. for S > 1, 
quantum fluctuations of the spin become less important. 



In the limit of large spin S, and at low excitation energies, it is natural to describe the 
ordered phase in terms of small fluctuations of the spins around their expectation values 
(cf. the description of the ordered phase of a crystal in terms of small fluctuations of the 
atoms around the ordered lattice sites). These fluctuations are conveniently represented 
in terms of spin raising and lowering operators: with = S^^iS^, it is straightforward 
to verify that 

K)^n] = ^mnS m , [<S+, S n ] = 2<5 mrl S'^. (2.15) 

Application of lowers (raises) the z-component of the spin at site m by one. To 

actually make use of the fact that deviations around |Q) are small, a representation 
known as the Holstein— Primakoff transformation 8 was introduced in which the spin 
operators S ,± ,5 are specified in terms of bosonic creation and annihilation operators a) 
and a: 



S m — a L(25' a ln a m) ^ ' S m — (2S aj^Qyjj) / a m , S m — S o^a. 



8 T. Holstein and H. Primakoff, Field dependence of the intrinsic domain magnetisation of a ferromag- 
net, Phys. Rev. 58, 1098 (1940). 
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Figure 2.8: Measurements of 
the spin-wave dispersion re- 
lations for the ferromagnet 
La .7Sr .3MnO 3 . 
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The utility of this representation is clear: When the spin is large S S> 1, an expansion 
in powers of 1/S gives S z m = S-ala m , S~ = (2S)^al+0(S- 1 / 2 ), and 5+ = (25)^ + 
0(S 1 / 2 ). In this approximation, the one-dimensional Hamiltonian takes the form 



H = 



- JiV5 2 + JSy~] jaj„a m + a^ +1 a m+ i - {a ] m a m+1 + h.c.) | + 0(5° 

m 

-JNS 2 + Sj>™+i - al)(a m+l - a m ) + O(S ). 



Bilinear in Bose operators, the approximate Hamiltonian can be diagonalised by Fourier 
transformation. With periodic boundary conditions, S m+ N = S m , a m+ N = a m , defining 



N 



a t = 



i 



B.Z. 



m— 1 



the Hamiltonian for the one dimensional lattice system takes the form (exercise) 



B.Z. 



H = 



-JNS 2 + ^La k ala k + O(S ) 



(2.16) 



where uj^ = 2JS{\ — cosfc) = 4 JS sin 2 (fc/2) represents the dispersion relation of the spin 
excitations. In particular, in the limit k — > 0, the energy of the elementary excitations 
vanishes, UJk — > JSk 2 (cf. Fig. 2.8). These massless low-energy excitations, known as 
magnons, describe the elementary spin-wave excitations of the ferromagnet. Taking 
into account terms at higher order in the parameter 1/5, one finds interactions between 
the magnons. 
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Quantum Antiferromagnet 

Having explored the low-energy excitation spectrum of the ferromagnet, we turn now to 
the spin S Heisenberg antiferromagnetic chain, 



H — J ^ S m ■ S m . j 



where J > 0. Such local moment antiferromagnetic phases frequently occur in the arena 
of strongly correlated electron systems. Although the Hamiltonian differs from its fer- 
romagnetic relative 'only' by a change of sign, the differences in the physics are dras- 
tic. Firstly, the phenomenology displayed by the antiferromagnetic Hamiltonian H de- 
pend sensitively on the geometry of the underlying lattice: For a bipartite lattice, i.e. 
one in which the neighbours of one sublattice A belong to the other sublattice B (cf. 
Fig. 2.9a), the ground states of the Heisenberg antiferromagnet are close 9 to a staggered 
spin configuration, known Neel 

Louis Neel 1904-2000: 1970 Nobel 
Laureate in physics for fundamental 
work and discoveries concerning anti- ' ,v 'v 
ferromagnetism and ferrimagnetism 
which have led to important appli- 
cations in solid state physics 



state, where all neighbouring spins 
are antiparallel. Again the ground 
state is degenerate, i.e. a global ro- 
tation of all spins by the same amount 
does not change the energy. By con- 
trast, on non-bipartite lattices such 




as the triangular lattice shown in 

Fig. 2.9b, no spin arrangement can be found wherein which each and every bond can 
recover the full exchange energy J. Spin models of this kind are said to be frustrated. 



O: A 
X : B 



O X O X O X O 

x O X O X O x 

O x O X O x O 

x o x o x o x 

O x o X O x o 
a) 



b) 



Figure 2.9: (a) Example of a two-dimensional bipartite lattice and (b) a non-bipartite lattice. 
Notice that, with the latter, no antiferromagnetic arrangement of the spins can be made that 
recovers the maximum exchange energy from each and every bond. 



> EXERCISE. Using only symmetry arguments, specify one of the possible ground states 
of a classical three site triangular lattice antiferromagnet. (Note that the invariance of the 
Hamiltonian under a global rotation of the spins means that the there is manifold of continuous 

9 It is straightforward to verify that the classical ground state — the Neel state — is now not an 
exact eigenstate of the quantum Hamiltonian. The true ground state exhibits zero-point fluctuations 
reminiscent of the quantum harmonic oscillator or atomic chain. However, when S 2> 1, it serves as a 
useful reference state from which fluctuations can be examined. 
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degeneracy in the ground state.) Using this result, construct one of the classical ground states 
of the infinite triangular lattice. 



Returning to the one-dimensional system, we first note that a chain is trivially bi- 
partite. As before, our strategy will be to expand the Hamiltonian in terms of bosonic 
operators. However, before doing so, it is convenient to apply a canonical transformation 
to the Hamiltonian in which the spins on one sublattice, say B, are rotated through 180° 
about the S B H> S%, S V B H> —S B , and S Z B )->■ — S Z B , i.e. when represented 

in terms of the new operators, the Neel ground state looks like a ferromagnetic state, 
with all spins aligned. We expect that a gradual distortion of this state will produce the 
antiferromagnetic analogue of the spin-waves discussed in the previous section. 

Represented in terms of the transformed operators, the Hamiltonian takes the form 



qz qz 
D m D m+l 



q+ q+ I q- q 

kJ m'~'m+l ' ^m^; 



m+1 J 



Once again, applying an expansion of the Holstein-Primakoff representation, S m ~ (2S) ' a' m , 
etc., one obtains the Hamiltonian 



H = -NJS 2 + 



t „t 



a m a m + a m+l a m+l + a m a m+l + a m a m+l 



O(S ) . 



At first sight the structure of this Hamiltonian, albeit bilinear in the Bose operators, looks 
akward. However, after Fourier transformation, a m = N~ l / 2 ^ k e~ tkm a k , it assumes the 
more accessible form (exercise) 



H = -NJS(S + 1) + JSJ2(4 a_ fc )( 7 1 fc 



0(S ), 



where 7 fc = cos k. 

Quadratic in the bosonic operators, the Hamiltonian can be again diagonalised by 
canonical transformation, i.e. 

N. N. Bogoliubov 1909-1992: Theoretical 
physcists aclaimed for his works in nonlinear 
mechanics, statistical physics, theory of super- 
fluidity and superconductivity, quantum field 
theory, renormalization group theory, proof of 
dispersion relations, and elementary particle 
theory. 



a transformation of the field 
operators that preserves the 
commutation relations. In 
the present case, this 



is 



achieved by a Bogoliubov 
transformation . 




cosh 6k 
- sinh 6k 



— sinh 6k 
cosh 6 k 



(2.17) 



> EXERCISE. Construct the inverse transformation. Considering the commutation relations 
of the operators a a , where a\ = a and 02 = a', explain why the Bogoliubov transformation is of 
the form of a Lorentz transformation. If operators a obeyed fermionic commutation relations, 
what form would the transformation take? 



Quantum Condensed Matter Field Theory 



2.2. APPLICATIONS OF SECOND QUANTISATION 



37 




(3/4,1/4) (1/2,1/2) (1/2,0) (3/4,1/4) (1,0) 



Figure 2.10: Experimentally obtained 
spin-wave dispersion of the high-T c par- 
ent compound LaCu04 — a prominent 
spin 1/2 antiferromagnet. Figure repro- 
duced from R. Coldea et ai, Phys. Rev. 
mm Lett. 86, 5377 (2001). 



Applying the Bogoliubov transformation, and setting tanh 29^ 
assumes the diagonal form (exercise) 



-7fc, the Hamiltonian 



H 



-NJS- 



sink] a\oi k 



(2.18) 



Thus, in contrast to the ferromagnet, the spin-wave excitations of the antiferromagnet 
exhibit a linear dispersion in the limit k — > 0. Surprisingly, although developed in the 
limit of large spin, experiment shows that even for 5=1/2 spin chains, the integrity of 
the linear dispersion is maintained (see Fig. 2.10). 



2.2.5 Bogoliubov theory of the weakly interacting Bose gas 

Ealier, we explored the influence of interactions on the electron gas. When interactions 
are weak, it was noted that the elementary collective excitations are reminiscent of the 
excitations of the free electron gas — the Fermi-liquid phase. When the interactions are 
strong, we discussed a scenario in which the electron liquid can condense into a solid 
insulating phase — the Mott transition. In the following section, we will discuss the 
properties of a quantum liquid comprised of Bose particles — the weakly interacting 
dilute Bose gas. 

Let us consider a system of N Bose particles confined to a volume L d and subject to 
the Hamiltonian 



H = 5^e k o£a k + i I d d xd d x' a t (x)a t (x / )F(x - x')a(x')a(x) 



where = and V(x) denotes a weak repulsive pairwise interaction. In the case 
of a Bose gas, this assumption is connected with the fact that, even for infinitesimal 
attractive forces, a Bose gas cannot stay dilute at low temperatures. In Fourier space, the 
corresponding two-body interaction can be expressed as (exercise) 

k,k',q 
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where ctk = J d d xe lkx a(x), and V q = J d d xe~* q ' x V(x). In the following, we will be 
interested in the ground state and low-lying excitations of the dilute system. In this case, 
we may distill the relevant components of the interaction and considerably simplify the 
model. 

In the ground state, the particles of an ideal (i.e. non-interacting) Bose gas condense 
into the lowest energy level. In a dilute gas, because of the weakness of the interactions, 
the ground state will differ only slightly from the ground state of the ideal gas, i.e. the 
number of particles iVo in the condensate will still greatly exceed the number of particles 
in other levels, so that JV - JV < JV. Since the number of particles in the condensate is 
specified by the number operator N = a\ =0 ak=o = O(N) S> 1, matrix elements of the 
Bose operators scale as a ~ O(^/N ). w This means that, from the whole sum in the 
interaction, it is sufficient to retain only those terms which involve interaction with the 
condensate itself. Taking V q = V constant, one obtains (exercise) 



a k a k + al k a_ k + - (a_ k a k + a k a - k 



+ 0(N? 

k/0 L 



. 



Terms involving the excited states of the ideal gas have the following physical inter- 
pretation: 

> VcL^dk represents the 'Hartree-Fock energy' of excited particles interacting with the 
condensate; 11 

> y(a_kflk + a k a -k) represents creation or annihilation of excited particles from the 
condensate. Note that, in the present approximation, the total number of particles 
is not conserved. 

Now, using the identity N = N + X^k/o a k a k to trade for N , the total Hamiltonian 
takes the form 

H = + ( e k 0) + Vn ) (4«k + a ! _ k a-k) + ^ (a-k«k + 4«-k) 

k/0 L 

where n = N / L d represents the total number density. This result may be compared 
with that obtained for the Hamiltonian of the quantum antiferromagnet in the spin 
wave approximation. Applying the Bogoluibov transformation (2.17): a k = cosh^k^k — 
sinh^ka^k' e ^ c -) with (exercise) 

smh 2 k =-(^- ij, 

where e k = [(e k °^ + Vn) 2 — (Vn) 2 ] 1 / 2 , one obtains 

k/0 k/0 



10 Note that the commutator [ao, aj] = 1 is small as compared to ao and a\ allowing the field operators 
to be replaced by the ordinary c-number vNq. 

11 Note that the contact nature of the interaction disguises the presence of the direct and exchange 
contributions. 
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From this result, we find that the spectrum of low energy excitations scales linearly as 
6k — ftc|k| where the velocity is given by c = (Vn/m) 1 / 2 while, at high energies (when 
k ~ ha = mc/h), the spectrum becomes free particle-like. 12 

> INFO. Since the number operator a k Ok can assume only positive values, one can infer 
the ground state wavefunction from the condition Ok|g- s -) — 0. Noting that the Bogoliubov 
transformation can be written as Ok = UajJJ -1 , with U = exp[^ k ^ ^(a k a -k — a kO—k)] 
(exercise), one can obtain the ground state as |g.s.) = U\&), where |<E>) = (a k =o) iV l^) denotes 
the ground state of the ideal Bose gas and | Q) the vacuum. The proof follows as 



= a k#0 |$> = U~ l UakU- 1 U\$). 

For the contact interaction, the corresponding ground state energy diverges and must be 'regu- 
larised'. 13 In doing so, one obtains Eq = — \ ^wo( e k°^ + n ^ — e k — ) which, when 
summed over k, translates to the energy density 



1? ~ ~T~ 



l + -l^(na 3 ) 1 / 2 

15V7T 



where a ~ {m/^-Kh 2 )V denotes the scattering length of the interaction. 

Finally, one may estimate the depletion of the condensate due to interaction. 



^^ = -^(g.s.|a k a k |g.s.) = -^smh k = - / sinh k = ^-k , 



k/0 k/0 



i.e. ca. one particle per "coherence length" £ ~ 1/fco- Recast using the scattering length, one 
obtains 

TV 3y^ 1 ' • 



How do these prediction compare with experiment? 14 When cooled to temperatures 
below AK, 4 He condenses from a gas into a liquid. The 4 He atoms obey Bose statistics 
and, on cooling still further, the liquid undergoes a transition to a superfluid phase in 
which a fraction of the Helium atoms undergo Bose-Einstein condensation. Within this 
phase, neutron scattering can be used to probe the elementary excitations of the system. 
Fig. 2.11 shows the excitation spectrum below the transition temperature. As predicted 
by the Bogoluibov theory, the spectrum of low-energy excitations is linear. The data also 
show the limitations of the weakly interacting theory. In Helium, the steric interactions 
are strong. At higher energy scales an important second branch of excitations known as 
rotons appear. The latter lie outside the simple hydrodynamic scheme described above. 

12 Physically, the effect of the interaction is to displace particles from the condensate even at T = 0. 
13 For details see, e.g., Ref. [1]. 

14 For a review of the history of Bose-Einstein condensation see, e.g. Griffin, cond-mat/9901123 
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Figure 2.11: The dispersion curve for 
4 He for the elementary excitations at 
1.1K (from Cowley and Woods 1971). 
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While studies of 4 He allow only for the indirect manifestations of Bose-Einstein con- 
densation (viz. superfluidity and elementary excitations of the condensate), recent investi- 
gations of dilute atomic gases allow momentum distributions to be explored directly. Fol- 
lowing a remarkable sequence of technological breakthroughs in the 90s, dilute vapours of 
alkali atoms, confined in magnetic traps, were cooled down to extremely low temperatures, 
of the order of fractions of microkelvins! Here the atoms in the vapour behave as quantum 
particles obeying Bose or 



Fermi statistics depend- 
ing on the atomic num- 
ber. By abruptly re- 
moving the trap, time-of- 
flight measurements al- 
low the momentum dis- 
tribution to be inferred 
directly (see Fig. 2.12a). 



Steven Chu 1948-, Claude 
Cohen- Tannoudji 1933- 
and William D. Phillips 

1948-: 1997 Nobel Laureates 
in Physics for development 
of methods to cool and trap 
atoms with laser light. 




Below a certain critical temperature, these measurements re- 
vealed the development of a sharp peak at low momenta in Bose gases providing a clear 
signature of Bose-Einstein condensation. 

In the condensed phase, one may measure the sound wave velocity by effecting a 
density fluctuation using the optical dipole force created by a focused, blue-detuned laser 
beam. By measuring the speed of propagation of the density fluctuation, the sound wave 
velocity can be inferred. The latter, shown in Fig. 2.12b, shows a good agreement with 
the theoretical prediction of the Bogoliubov theory. 



2.3 Summary 

This concludes our discussion of the second quantisation and its applications to problems 
in many-body quantum mechanics. Beyond qualitative discussions, the list of applications 
encountered in this chapter involved problems that were either non-interacting from the 
outset, or could be reduced to a quadratic operator structure by a number of suitable 
manipulations. However, we carefully avoided dealing with interacting problems where 
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Figure 2.12: Images of the velocity distribution of rubidium atoms by Anderson et al. (1995), 
taken by means of the expansion method. The left frame corresponds to a gas at a temperature 
just above condensation; the center frame, just after the appearance of the condensate; the right 
frame, after further evaporation leaves a sample of nearly pure condensate. The field of view is 
200/im x 200/im, and corresponds to the distance the atoms have moved in about l/20s. Speed 

of sound, c vs. condensate peak density Nq for waves propagating along the axial direction in the 

l li 

condensate. Data taken from Kurn et al. (1997) compared to theoretical prediction c ~ N ' . 

no such reductions are possible. Yet it should be clear already at this stage of our dis- 
cussion that completely or nearly solvable systems represent only a small minority of the 
systems encountered in condensed matter physics. What can be done in situations where 
interactions, i.e. operator contributions of fourth or higher order, are present? 

Generically, interacting problems of many-body physics are either fundamentally inac- 
cessible to perturbation theory, or they necessitate perturbative analyses of infinite order 
in the interaction contribution. Situations where a satisfactory result can be obtained by 
first or second order perturbation theory are exceptional. Within second quantisation, 
large order perturbative expansions in interaction operators leads to complex polynomials 
of creation and annihilation operators. Quantum expectation values taken over such struc- 
tures can be computed by a reductive algorithm, known as Wick's theorem. However, 
from a modern perspective, the formulation of perturbation theory in this way is not very 
efficient. More importantly, problems that are principally non-perturbative have emerged 
as a focus of interest. To understand the language of modern quantum condensed matter, 
we thus need to develop another layer of theory, known as field integration. However, 
before discussing quantum field theory, we should understand how the concept works in 
principle, i.e. on the level of point particle quantum mechanics. This will be the subject 
of the next chapter. 
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